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FINAL  REPORT:  Contract  W91  INF- 11-1-0481,  “Not  different,  Just  Better:  the  adaptive 
evolution  of  an  enzyme” 

1)  Summary 

We  know  much  about  the  consequences  of  adaptation,  but  little  of  the  underlying  molecular 
causes.  Our  program  combined  biochemical,  biophysical  and  evolutionary  experiments  to 
rigorously  examine  the  molecular  and  biochemical  basis  of  adaptation,  mediated  by  changes 
in  pyruvate  kinase  found  in  Richard  Lenski’s  E.  coli  long-term  evolution  experiment.  This 
program  was  aimed  at  uncovering  the  molecular  basis  for  a  series  of  adaptive  mutations  in  a 
key  allosteric  enzyme.  We  chose  the  enzyme  pyruvate  kinase  (gene  pykF )  from  Escherichia 
coli  as  our  model  system,  since  1 )  it  catalyzes  a  reaction  central  to  energy  metabolism,  2)  it  is 
precisely  regulated  by  allostery  and  the  adaptation  of  allostery  is  unknown,  and  3)  multiple 
experiments  by  others  have  demonstrated  that  adaptive  mutations  reproducibly  occur  in  this 
enzyme  when  E.  coli  is  grown  under  low  glucose  conditions. 

We  have  demonstrated,  for  the  first  time,  that  all  the  pykF  mutations  found  in  the  E.  coli  long¬ 
term  evolution  experiment  confer  a  significant  adaptive  fitness  advantage  to  the  bacterium. 
Moreover,  we  have  shown  that  the  mutations  are  likely  to  have  a  similar  proximal  basis  of  the 
benefit  they  confer:  a  shorter  lag  phase  of  the  growth  cycle.  Not  only  are  adaptive  mutations 
often  parallel  at  the  physiological  and  genetic  level,  they  are  also  parallel  at  the  molecular 
level.  By  measuring  the  effect  of  pykF  mutations  in  a  range  of  resources  that  have  different 
characterized  uptake  mechanisms,  we  have  made  strides  in  determining  the  underlying 
molecular  basis  of  their  benefit. 

The  key  finding  from  our  program  is  that  while  the  adaptive  mutations  in  pyruvate  kinase 
significantly  alter  the  catalytic  and  allosteric  function  of  the  enzyme,  they  do  not  change  the 
overall  structure  of  the  proteins  or  their  thermal  stability.  Instead,  we  have  demonstrated 
(using  cutting  edge  deuterium  exchange  studies  and  molecular  dynamic  simulations)  that  the 
dynamics  of  key  loops  surrounding  the  active  site  are  altered,  which  explain  the  changes  in 
catalytic  function.  Thus,  we  propose  that  protein  dynamics  is  an  important  and  general 
underlying  molecular  cause  of  adaptation  in  enzymes. 

Serendipitously,  we  solved  the  structure  of  E.  coli  pyruvate  kinase  A301S  with  and  without  its 
allosteric  effector  (fructose-1 ,6-bisphosphate),  allowing  us  to  propose  a  mechanism  of 
allosteric  regulation  in  this  enzyme.  In  particular,  the  binding  of  fructose-1, 6-bisphosphate 
causes  a  change  in  the  relative  orientation  of  domains  within  the  monomer,  which  in  turn 
affects  the  dynamics  of  key  active  site  loops  involved  in  substrate  binding.  Deuterium 
exchange  studies  support  our  hypothesized  mechanism,  and  also  demonstrate  that  the 
allosteric  signal  is  communicated  from  the  allosteric  binding  site  to  the  active  site  by 
destabilization  of  the  beta-strands  that  form  the  barrel  in  the  (a/(3)e-barrel  active  site  domain. 


2)  We  have  achieved  each  of  our  stated  goals 

Question  1:  How  do  the  effects  of  adaptive  mutations  depend  on  their  genetic  background? 

•  Aim  1)  Allelic  replacement  of  the  evolved  and  ancestral  alleles  in  the  ancestral  and 
evolved  genetic  backgrounds.  Achieved 

•  Aim  2)  Quantifying  and  comparing  the  fitness  of  the  various  pykF  alleles  in  different 
genetic  backgrounds.  Achieved 
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Question  2:  What  is  the  repeatability  of  molecular  evolution?  That  is,  are  adaptive  mutations 
in  the  same  gene,  but  replicate  populations,  functionally  parallel? 

•  Aim  3)  Expression,  purification  and  functional  analysis  of  evolved  pyruvate  kinase  type 
1  enzymes.  Achieved 

•  Aim  4)  Investigating  the  structural  changes  in  the  evolved  pyruvate  kinase  type  1 
enzymes.  Achieved 

•  Aim  5)  Unraveling  the  metabolic  consequences  of  the  evolved  pyruvate  kinase  type  1 
enzymes.  Achieved 

Question  3:  Do  the  adaptive  pyruvate  kinase  mutations  have  secondary  effects? 

•  Aim  6)  Determine  the  fitness  of  evolved  pykF  alleles  in  novel  resources.  Achieved 

Thesis  and  Manuscript  preparation  -  Driven  by  PI  Dobson  and  Cooper.  On  Track 


The  ARO  contract  has  supported  the  following  staff  and  students: 

•  1  Postdoctoral  Fellow  (Dr  Sarah  Kessans,  University  of  Canterbury) 

•  2  PhD  students  (now  Dr  Fen  Peng,  University  of  Houston;  Dr  Katherine  Donovan,  UC) 

•  3  Honours  students  (Elena  Sugrue,  UC;  Rosie  Shaw,  UC;  Kristina  Duan,  UoH) 

•  1  Summer  studentship  (Anji  Stampfli,  UC). 
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D.W.  Griffin,  Geoffrey  B.  Jameson  and  Renwick  C.J.  Dobson  (2015)  Grappling  with  anisotropic 
data,  pseudo-merohedral  twinning  and  pseudo-translational  non-crystallographic  symmetry: 
A  case  study  involving  pyruvate  kinase.  Provisionally  accepted  Acta  Crystallographies  D: 
Biological  Crystallography 

2)  Katherine  A.  Donovan,  Shaolong  Zhu,  Peter  Liuni,  Fen  Peng,  Sarah  Kessans,  Tim  F.  Cooper,  Derek 
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Dancing  to  different  tune:  Natural  selection  of  an  enzyme  by  adaptation  fine-tunes  protein 
dynamics.  To  be  drafted  and  will  be  submitted  to  Science. 
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The  dynamic  lid  domain  of  pyruvate  kinase:  not  involved  in  allostery,  but  critical  for  reaction 
promiscuity.  To  be  drafted  based  on  Elena  Sugrue’s  Honours  report  and  submitted  to  FEBS 
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Introduction 

How  do  organisms  adapt  to  novel  environments?  Adaptation  is  the  process  by  which  a 
population  moves  towards  a  phenotype  that  represents  a  better  fit  to  the  environment.  We 
know  a  great  deal  about  the  consequences  of  adaptation;  for  example,  the  speed  with  which 
adaptive  mutations  of  certain  fitness  effects  will  spread  in  a  population  and  how  this  will  be 
affected  by  changing  population  parameters  (size,  recombination,  etc).  In  contrast,  we  know 
very  little  about  the  underlying  molecular  causes  of  adaptation;  that  is,  how  mutations  act  to 
produce  an  organism’s  phenotype.  Indeed,  with  few  exceptions,  such  as  the  emergence  of 
drug  resistance,  the  link  between  genetic  changes  and  their  effect  on  fitness  has  been  left  as 
a  black  box — mutations  go  in,  the  phenotype  comes  out.  Answers  to  fundamental  questions 
are  not  just  unknown,  but  unexplored.  Understanding  the  mapping  between  the  molecular 
consequences  of  adaptive  mutations  and  their  phenotypic  consequences  will  allow  us  to  ask 
new  types  of  questions:  1 )  How  do  the  effects  of  an  adaptive  mutation  depend  on  genetic 
background?  2)  Are  independent  adaptive  mutations  in  the  same  gene  functionally  parallel; 
that  is,  do  they  confer  the  same  function  to  the  evolved  protein?  3)  What  role  does  protein 
stability  play  in  constraining  the  ability  of  an  enzyme  to  adapt? 

Testing  models  of  adaptation  poses  several  challenges — not  least  the  long  time  scales  that 
can  be  required  for  adaptation  to  occur.  To  circumvent  this  limitation  we  build  on  a  long-term 
laboratory  study  in  which  12  independent  replicate  E.  coli  populations  have  been  evolved  in  a 
glucose-limiting  environment  for  60,000  generations  (~1  million  human  years!).  During  this 
time,  the  replicate  populations  increased  in  fitness  relative  to  their  common  ancestor. 
Remarkably,  the  genetic  changes  underlying  this  adaptation  were  often  parallel,  being 
concentrated  in  relatively  few  genes.  In  particular,  mutations  in  one  gene,  pykF,  occurred 
independently  (and  were  fixed)  in  all  12  long-term  populations — a  signature  that  they  are 
likely  to  be  adaptive.  Moreover,  in  one  case  the  same  mutation  in  pykF  occurred 
independently  in  three  populations.  Our  research  program  studied  why  pykF  serves  as  an 
adaptive  focus.  To  do  this,  we  have  developed  a  molecular  ‘picture’  of  the  effect  of  the 
evolved  changes  on  the  structure  and  function  of  the  gene  product — pyruvate  kinase. 

Pyruvate  kinase  catalyzes  a  key  regulatory  step  in  glycolysis — the  biochemical  pathway 
central  to  energy  metabolism.  The  catalyzed  reaction  involves  transfer  of  a  phosphoryl  group 
from  phosphoenolpyruvate  to  ADP,  generating  the  important  biomolecules  pyruvate  and  ATP. 
Tight  regulation  of  glycolysis  is  crucial,  since  gluconeogenesis  (the  opposing  process) 
would  otherwise  occur  simultaneously,  leading  to  wasteful  cycling  of  biochemical 
intermediates.  We  hypothesized  that  pyruvate  kinase  is  a  focal  point  for  parallel  evolution 
because  it  is  a  key  control  point  in  energy  metabolism.  Although  the  adaptive  mutations  in 
PK1  are  distributed  throughout  the  gene  sequence,  when  mapped  onto  the  pyruvate  kinase 
structure  they  clustered  to  three  hot  spots.  We  predicted  that  the  mutations  would  have 
substantial  effects  on  the  catalytic  activity  and/or  allosteric  regulation  of  the  evolved  enzymes, 
since  the  currently  accepted  mechanism  of  PK1  regulation  occurs  via  rearrangement  of  the 
tetramer  interfaces. 

Understanding  and  predicting  the  effect  of  mutations  at  the  organismal  level  is  a  major 
challenge  of  evolutionary  biology.  Empirical  data  addressing  why  some  mutations  are 
adaptive  (and,  conversely,  why  most  are  not)  will  have  major  implications  for  our  ability  to 
predict  and  understand  the  outcome  of  evolution,  both  in  vitro  and  in  vivo. 
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Results  to  proposed  questions  and  data. 

Question  1:  How  do  the  effects  of  an  adaptive  mutation  depend  on  genetic 
background? 

Summary: 

We  have  discovered  that  the  magnitude  of  the  fitness  benefit  conferred  by  the  adaptive 
mutations  in  pyruvate  kinase  depends  on  the  background  on  which  they  arise,  although  they 
are  almost  always  positive.  This  strongly  suggests  that  pyruvate  kinase  is  a  generalist  (in  the 
sense  that  evolved  mutations  can  confer  benefits  in  many  backgrounds),  but  is  also  involved 
in  epistatic  interactions.  Moreover,  we  discovered  that  the  different  adaptive  mutations  are 
parallel  in  their  effect,  providing  the  first  evidence  that  although  the  mutations  are  spread 
across  the  gene,  their  fitness  effect  is  very  likely  to  be  the  same. 

For  extended  details  see  Fen  Peng’s  PhD  thesis,  Chapter  2. 

•  Aim  1)  Allelic  replacement  of  the  evolved  alleles  in  the  ancestral  genetic  background. 

Achieved 

Detailed  Results: 

We  constructed  two  series  of  strains  that  were  designed  to  allow  comparison  of  the  effect  of 
independently  evolved  mutations  in  pykF  in  a  common  genetic  background  (the  ancestor  to 
our  evolution  experiment)  and  in  the  evolved  background  in  which  they  are  present  in  the 
evolution  experiment.  We  found  that  the  different  point  mutation  alleles  of  pykF  had  similar 
fitness  effects  in  the  ancestral  background  (Figure  1).  Interestingly,  these  effects  were 
significantly  different  from  that  of  a  deletion  allele,  clearly  indicating  that  the  mutant  enzymes 
retain  some  function.  The  fitness  of  the  evolved  pykF  genes  added  to  the  ancestor  increased 
in  the  range  from  6%  to  12%  (see  Figure  1 ). 

Effects  were  generally  higher  and  more  variable  in  the  evolved  backgrounds,  where  they 
ranged  from  0  to  25%  (Figure  2).  We  draw  one  important  conclusions  from  these  results.  The 
effects  of  focal  mutations  are  not  constant,  but  depend  on  the  presence  of  additional 
mutations  as  they  arise  during  evolution. 

•  Aim  2)  Quantifying  and  comparing  the  fitness  of  the  various  pykF  alleles  in  different 
genetic  backgrounds.  Achieved 

Detailed  Results: 

We  tested  two  explanations  for  our  finding  that  the  evolved  pykF  enzymes  had  different 
fitness  effects  in  the  ancestor  and  in  their  evolved  background  (Figure  3):  1 )  they  are  caused 
by  different  pyruvate  kinase  enzymes,  or  2)  they  are  caused  by  the  different  genetic 
backgrounds  in  which  they  occur.  To  do  this,  we  moved  each  of  three  pykF  alleles  (WT, 
A301S,  Del)  into  1 1  different  evolved  backgrounds  (Figure  4).  We  found  that  each  of  the 
transferred  alleles  had  similar  effects  in  the  same  background  and  very  different  effects 
across  backgrounds. 

Additional  work: 

The  evolved  alleles  result  in  a  shorter  lag  phase  growth.  As  a  first  step  in  determining  the 
mechanism  underlying  fitness  effects  of  mutations  in  pykF,  we  have  analyzed  growth  curves 
of  the  series  of  strains  made  by  adding  the  mutations  in  to  the  common  ancestral  background. 
We  find  a  significant  trend  toward  a  decrease  in  lag  time  -  the  time  taken  for  cells  to  begin 
growth  following  transfer  to  fresh  medium  (Figure  3A)  -  and  a  smaller  and  non-significant 
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trend  toward  faster  growth  rates  (Figure  3B).  These  results  are  consistent  with  previous  work 
proposing  that  pykF  mutations  might  provide  a  benefit  to  cells  by  allowing  them  to  maintain 
higher  PEP  (phosphoenolpyruvate)  -  the  substrate  of  the  pyruvate  kinase  -  levels  during 
stationary  phase,  increasing  the  rate  at  which  glucose  can  be  imported  into  cells  when  it  is 
added  to  the  environment. 
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Figure  1:  Fitness  of  ancestral  and  evolved 
pykF  genes  in  the  ancestral  genetic 
background.  Symbols  indicate  mean  and 
error  bars  indicate  the  95%  Cl  of  three 
replicate  fitness  measurements. 
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Figure  2:  Fitness  of  evolved  pykF  genes 
in  the  evolved  backgrounds.  Symbols 
indicate  mean  and  error  bars  indicate  the 
95%  Cl  of  three  replicate  fitness 
measurements. 


Genotype 
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Figure  3:  Lag  time  (A.)  and  growth  rates  (B.) 
of  the  evolved  pykF  allele  in  the  ancestor 
strain,  compared  directly  with  the  ancestor. 

Symbols  indicate  mean  and  error  bars  indicate 
the  95%  Cl  of  three  replicate  fitness 
measurements.  There  is  a  significantly  shorter 
time  for  evolved  strains,  except  for  the  deletion. 


Figure  4.  Fitness  effects  of  different 
pyruvate  kinase  enzymes  in  11 
independently  evolved  genetic 
backgrounds.  Different  symbol  colors 
indicate  the  effect  of  each  of  three  pykF 
genes  in  each  genetic  background.  Genetic 
backgrounds  are  clones  isolated  from 
populations  started  from  a  common  ancestor 
and  evolved  independently  for  20,000 
generations.  Symbols  indicate  mean  and 
error  bars  indicate  the  95%  Cl  of  three 
replicate  fitness  measurements. 
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Question  2:  Are  adaptive  mutations  in  the  same  gene,  but  replicate  populations, 
functionally  parallel? 

Summary: 

Following  Question  1,  in  which  we  discovered  that  the  adaptive  mutations  confer  a  parallel 
fitness  advantage,  we  asked  whether  the  evolved  enzymes  where  themselves  functionally 
altered  in  parallel  (i.e.  catalysis  and  regulation).  We  discovered  that  while  the  catalytic  and 
allosteric  function  of  the  evolved  enzymes  was  surprisingly  different,  the  overall  structure  was 
surprisingly  unaltered.  We  then  hypothesized  that  the  dynamics  of  the  enzymes  would  be 
different.  We  tested  our  hypothesis  using  hydrogen/deuterium  exchange  experiments,  which 
demonstrates  that  the  dynamics  of  the  proteins  are  indeed  different.  Thus,  we  have 
demonstrated  that,  in  this  case,  evolution  alters  enzyme  dynamics.  Lastly,  consistent  with  the 
central  role  of  pyruvate  kinase  in  energy  metabolism,  we  have  shown  that  the  metabolome  is 
significantly  altered  by  the  adaptive  mutations  in  pyruvate  kinase.  In  particular,  we  have 
discovered  that  the  TCA  cycle  and  fatty  acid  biosynthesis  are  down  regulated.  We  propose 
that  the  mutations  optimize  the  metabolome  for  growth  in  a  low  glucose  environment. 

For  extended  details  see  Katherine  Donovan’s  PhD  thesis,  Chapters  4,  5  and  6. 

•  Aim  2)  Expression,  purification  and  functional  analysis  of  evolved  pyruvate  kinase  type  1 
enzymes.  Achieved 

Detailed  Results: 

Figure  5  demonstrates  that  the  there  are  significant  differences  in  the  kinetic  properties  of  the 
eight  adaptive  mutations  when  compared  to  the  WT  enzyme  (i.e.,  the  wild  type  ancestral 
enzyme).  Importantly,  the  major  effect  is  on  PEP  binding  and  interestingly,  the  enzymes 
show  altered  allostery  with  respect  to  both  PEP  and  FBP.  This  is  an  exciting  result  and  may 
represent  the  first  observation  of  enzymes  altering  allostery  to  adapt  to  a  new  environment. 

Historically,  the  view  was  that  most  mutations  could  be  tolerated  without  loss  of  stability  or 
function.  A  recent  review  highlighted  a  change  in  this  view,  with  the  rationalization  that  most 
mutations  cause  some  change  in  stability,  usually  destabilizing.  It  has  been  proposed  that 
some  proteins  contain  inherent  robustness  (the  threshold  robustness  model),  such  that 
adaptive  mutations  can  be  tolerated  without  compromising  stability.  Given  that  the  adaptive 
mutations  accumulated  in  relatively  few  genes  in  the  LT  experiment,  and  in  particular  in 
pyruvate  kinase,  the  thermal  stability  of  wild  type  pyruvate  kinase  was  compared  to  the 
evolved  enzymes  (Figure  13).  These  studies  show  that  the  inherent  thermo-stability  of  the 
pyruvate  kinase  enzyme  is  rather  high  at  62  °C  and  that  the  adaptive  mutations  do  not  have  a 
significant  affect  on  thermal  stability.  Our  data  are  consistent  with  the  threshold  robustness 
model. 

•  Aim  3)  Investigation  of  structural  changes  of  the  evolved  pyruvate  kinase  type  1 
enzymes.  Achieved 

Detailed  Results: 

Given  that  several  of  the  mutations  (A301S,  A301T  and  G381A)  reside  at  the  tetrameric 
interface,  we  investigated  whether  these  adaptive  mutations  affected  the  oligomeric  state  of 
the  enzymes.  To  examine  the  association/dissociation  limit  of  wild  type  and  the  evolved 
pyruvate  kinase  enzymes,  analytical  ultracentrifugation  experiments  were  performed  (Figure 
6).  None  of  the  mutant  enzymes  had  altered  quaternary  function,  even  at  low  nM 
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concentrations.  These  data  demonstrate  that  the  adaptive  mutations  do  not  cause  changes 
in  the  quaternary  structure. 

To  examine  the  structural  changes,  we  initiated  atomic  resolution  crystallization  studies  on  the 
mutant  enzymes  (resolution  varying  from  2-2.6  A)  comparing  them  to  a  wild-type  structure 
that  we  determined  in  house.  Surprisingly,  each  of  the  evolved  enzymes  overlay  very  closely 
with  the  wild  type  structure  (Figure  7)  and  only  very  small  changes  could  be  found.  This  was 
surprising,  given  the  very  different  kinetic  and  allosteric  functions  of  the  enzymes.  We  verified 
that  in  solution  the  mutation  were  not  altering  the  structure  of  the  protein  using  small  angle  X- 
ray  scattering. 

•  Aim  4)  Unraveling  the  metabolic  consequences  of  the  evolved  pyruvate  kinase  type  1 
enzymes.  Achieved 

Finally,  we  collaborated  with  Dr  Silas  Villas-Boas  (University  of  Auckland)  to  explore  how 
changes  in  pyruvate  kinase  affect  the  metabolome.  The  most  significant  changes  in  the 
metabolome  are  shown  in  Figure  8.  Large  increases  in  the  concentrations  of  the  extracellular 
TCA  metabolites  (succinate,  malic  acid,  and  fumaric  acid)  and  fatty  acids,  suggest  that  flux 
through  this  pathway  is  low.  Interestingly,  although  the  changes  are  generally  consistent 
across  the  mutations,  there  are  significant  differences  with  two  mutations  (P70T  and  P70Q, 
Figure  8B),  which  suggest  that  there  is  more  than  one  way  to  optimize  the  metabolome. 
Overall,  we  interpreted  the  changes  in  the  metabolome  as  a  global  optimization  of 
metabolomic  function  to  the  low  glucose  environment.  It  is  unsurprising  that  pyruvate  kinase 
is  an  adaptive  focus  to  this  evolutionary  pressure,  given  its  central  role  in  metabolism. 

Additional  work: 

The  adaptive  mutations  alter  protein  dynamics.  Since  the  very  different  kinetics  of  the 
evolved  pyruvate  kinase  enzymes  (Figure  5),  cannot  be  explained  by  the  surprisingly  similar 
structures  when  compared  to  the  wild-type  enzyme,  we  hypothesized  that  the  changes  in 
kinetics  were  caused  by  the  mutations  altering  the  dynamics  of  the  proteins.  In  collaboration 
with  A/Prof  Derek  Wilson  (University  of  York,  Canada),  hydrogen/deuterium  exchange  mass 
spectrometry  studies  were  undertaken  to  gauge  whether  the  dynamics  of  the  evolved 
enzymes  was  different  to  the  wild-type  enzyme.  The  results  for  these  studies  are  consistent 
with  our  hypothesis:  global  deuterium  uptake  of  the  evolved  enzymes  is  different  to  each 
other  and  to  the  wild  type  enzyme  suggesting  that  they  have  different  global  dynamics  (Figure 
9).  Using  pepsin  digests,  we  have  shown  that  the  pattern  of  deuterium  exchange  across  the 
protein  is  also  different  (Figure  10).  Supporting  this  finding  are  molecular  dynamic  simulation 
studies  (Figure  1 1 ).  This  work  is  ongoing,  but  provides  strong  evidence  that  the  evolved 
mutations  are  altering  protein  dynamics. 

Unraveling  the  mechanism  of  allosteric  activation  of  E.  coli  pyruvate  kinase  by  FBP. 

Understanding  how  the  mutations  altered  the  allosteric  mechanism,  as  seen  in  the  kinetic 
studies  (Figure  5),  required  the  structure  of  pyruvate  kinase  with  bound  FBP.  Soaking  wild- 
type  and  A301T  pyruvate  kinase  crystals  with  the  allosteric  effector  caused  immediate 
cracking  of  the  crystals  and  no  diffraction  data  could  be  collected.  This  phenomenon  often 
occurs  when  the  ligand  binds  to  the  protein,  causing  significant  rearrangements  within  the 
protein,  such  that  the  protein  interactions  necessary  to  form  the  lattice  of  the  crystal  are 
broken.  In  contrast,  pyruvate  kinase  A301S  when  soaked  with  FBP  did  not  cause  crystal 
cracking  and  diffraction  could  be  collected  to  2.6  A.  This  structure  (Figure  12)  shows  that 
upon  binding  of  FBP,  the  orientation  of  the  monomers  within  the  tetramer  is  altered  with 
respect  to  one  another.  Each  twist  by  about  10  degrees  and  the  interfaces  of  the  tetramer  are 
significantly  changed.  When  we  compared  the  active  sites  of  the  FBP  bound  structure  with 
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the  unbound  structure,  we  notice  very  few  differences,  although  two  active  site  loops  become 
very  dynamic.  This  suggests  to  us  that  FBP  binding  remodels  the  active  site  to  better 
accommodate  PEP,  and  in  so  doing  activates  the  enzyme.  Importantly,  this  result  is 
consistent  with  our  finding  that  the  dynamics  of  the  active  site  loops  are  critical  for  function. 
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Figure  5:  Enzyme  kinetic  plots  of  the  wild  type  and  mutant  pyruvate  kinase  enzymes  demonstrate  that 
they  are  functionally  very  different  from  one  another.  These  plots  show  the  rate  of  reaction  relative  to 
increasing  phosphoenol  pyruvate  (PEP)  concentrations.  The  Ku,  kcaU  and  allosteric  effect  of  FBP  are  different 
across  this  set  of  enzymes. 
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Figure  6:  Analytical  ultracentrifugation  data  demonstrates 
that  the  pyruvate  kinase  A301T  is  a  tetramer  even  at  low 
nM  concentrations.  The  c(s)  distribution  (lower  panel)  shows 
only  a  single  peak  that  is  the  same  as  the  wild-type  enzyme. 
The  upper  panels  show  the  raw  analytical  ultracentrifugation 
data. 
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Figure  7:  The  crystal  structures  of  each  of 
the  evolved  pyruvate  kinase  enzymes 
compared  to  the  wild-type  enzyme  reveal 
that  the  evolved  enzymes  are  very  similar  in 
atomic  structure.  The  resolution  of  the  evolved 
pyruvate  kinases  ranged  from  2.0  -  2.6  A.  The 
Rfree  statistics  ranges  from  18%  -  25%.  The 
adaptive  mutation  could  be  resolved  in  the 
electron  density. 
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Figure  8.  Pyruvate  kinase  mutations  alter  the  cellular  metabolome.  A.  Significant  changes  in  the  relative 
abundance  of  extracellular  metabolites  were  focused  on  TCA  cycle  substrates  (malic  acid,  fumaric  acid,  2- 
oxoglutaric  acid,  and  succinic  acid;  illustrating  an  increase  in  relative  abundance  in  mutants  compared  to  wild 
type  cells),  fatty  acids  (linoleic  acid;  varying  across  mutants  and  wild  type),  and  putrescine  (varying  across 
mutants,  but  generally  lower  in  abundance  as  compared  to  the  wild  type).  B.  Significant  changes  in  the  relative 
abundance  of  intracellular  metabolites  were  clustered  around  fatty  acids  (generally  lower  abundances  in  mutants 
as  compared  to  wild  type),  TCA  cycle  substrates  (cis-aconitic  acid  higher  in  mutants,  succinic  acid  lower  in 
mutants),  and  amino  acids  (varied  across  mutants  and  wild  type). 
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Figure  9.  Global  hydrogen/deuterium 
exchange  studies  for  the  pyruvate  kinase 
enzymes  demonstrate  that  the  mutations 
cause  a  change  in  protein  dynamics.  This 
panel  plots  the  uptake  of  deuterium  versus 
time.  The  mutants  show  different  initial  uptake 
rates  and  different  total  uptake  (plateau). 
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Figure  10.  Peptide 
hydrogen/deuterium  exchange 
studies  for  the  pyruvate  kinase 
enzymes  demonstrate  that  the 
mutations  cause  a  change  in 
protein  dynamics.  The  colours 
indicate  the  percentage  deuterium 
uptake  after  2  s  (plateau  region  in 
Figure  9).  Note  that  the  pattern  of 
uptake  is  different  between  the  two 
proteins,  demonstrating  that  the 
dynamics  is  different. 


Figure  11.  Molecular  dynamic 
simulations  (500  ns)  of  the  wild  type 
and  A301T  pyruvate  kinase  enzymes 
demonstrate  differences  in  protein 
dynamics.  In  this  figure,  the  active  site 
domain  (a/|3)8-barrel  (A)  and  allosteric 
binding  domain  (B)  illustrate  broad 
changes  in  dynamics.  The  (a/|3)8-barrel 
of  the  wild  type  is  considerably  more 
stable  than  that  of  the  A301T  mutant, 
whereas  the  allosteric  binding  domain  is 
more  stable  in  the  wild-type  protein 
compared  to  the  mutant  enzyme. 
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Figure  12.  Overlay  of  the  A301S  pyruvate  kinase  with  FBP 
(blue)  and  without  FBP  (green)  showing  the  change  in 
conformation  that  accompanies  FBP  binding  (FBP  shown 
as  red  spheres).  From  these  data,  we  are  able  to  examine 
the  mechanism  of  FBP  binding  and  ‘see’  how  binding 
changes  the  conformation  of  the  active  site.  Consistent  with 
our  hydrogen/deuterium  exchange  studies,  loops  at  the  active 
site  (N-terminal  of  the  (a/(3)8-barrel)  have  higher  B-factors 
when  FBP  is  bound,  suggesting  that  they  are  more  dynamic. 
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Figure  13.  1-Anilino-8-naphthalene  sulfonate  (ANS)  binding  assays  demonstrate  that  the  evolved 
enzymes  are  approximately  equal  in  stability  when  compared  to  the  wild-type  enzyme  (blue).  While  there 
is  some  changes,  e.g.  G381A  (grey)  has  a  lower  melting  temperature  (54  °C)  as  compared  to  the  wild-type 
enzyme  (blue,  62  °C),  overall  the  mutations  do  not  consistently  alter  the  Tm  (i.e.  all  lower  or  all  higher). 
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Question  3:  Do  the  adaptive  pyruvate  kinase  mutations  have  secondary  effects? 

Summary: 

Even  when  the  fitness  effect  of  an  adaptive  mutation  can  be  estimated,  it  has  usually  not 
been  determined  how  the  mutation  causes  the  fitness  effect  that  it  does.  We  previously 
demonstrated  that  independently  selected  adaptive  mutations  in  pykF  conferred  similar 
fitness  effects  when  measured  in  a  defined  background.  It  is  unknown,  however,  if  these 
parallel  mutations  affect  fitness  through  the  same  underlying  mechanism.  Although  this  may 
seem  likely,  biochemical  studies  of  the  evolved  pyruvate  kinase  type  1  enzymes  reveal  a  wide 
array  of  changes  in  enzyme  kinetics  that  include  increases  and  decreases  in  activity  and 
stability  relative  to  the  wild  type  enzyme,  so  that  it  cannot  be  assumed  that  they  have  the 
same  mechanism  of  action.  To  examine  the  basis  of  the  benefit  conferred  by  mutations  in 
pyruvate  kinase  type  1,  we  proposed  and  tested  a  mechanistic  hypothesis  whereby  mutations 
in  pyruvate  kinase  act  to  increase  stationary  phase  concentrations  of  the  chemical  PEP, 
which  in  turn  facilitates  more  rapid  uptake  of  glucose  when  cells  are  transferred  to  fresh 
medium.  Consistent  with  predictions  of  this  hypothesis:  (1 )  pykF  mutations  tended  to  confer 
greater  benefits  in  substrates  requiring  PEP  for  uptake  than  in  alternative  substrates  (Figure 
14),  (2)  stationary  phase  PEP  levels  are  lower  in  the  wild  type  strain  than  in  strains  containing 
any  of  the  different  evolved  pyruvate  kinase  type  1  enzymes  (Figure  15),  and  (3)  the  pykF 
mutations  confer  most  of  their  benefit  through  shortening  the  lag  time  inherent  in  each  daily 
growth  cycle,  not  by  increasing  cell  growth  rate  (Figure  16).  Our  results  suggest  that,  while 
parallel  mutations  change  the  in  vitro  biochemical  properties  of  pyruvate  kinase  type  1  in  a 
range  of  different  ways,  the  net  effect  of  these  changes  is  to  confer  similar  adaptive  changes 
on  cell  metabolism. 

For  extended  details  see  Fen  Peng’s  PhD  thesis,  Chapters  3  and  4. 

•  Aim  6)  Determine  the  fitness  of  evolved  pykF  alleles  in  novel  resources.  Achieved 
Detailed  Results: 

PEP  is  not  only  the  substrate  of  pyruvate  kinase  type  1 ,  but  is  also  used  as  a  source  of 
energy  to  drive  the  phosphotransferase  system  (PTS),  a  system  required  for  the  transfer  of 
some  chemicals  into  the  cell.  If  mutations  in  pykF  affect  stationary  phase  levels  of  PEP,  we 
expect  them  to  have  a  general  effect  on  the  initial  translocation  of  PTS  resources  and 
therefore  for  pykF  mutants  to  generally  have  higher  fitness  in  PTS  than  non-PTS  resource 
environments.  We  note  that  pykF  mutations  may  well  have  pleiotropic  effects  that  cause 
fitness  to  vary  across  environments,  so  this  prediction  does  not  imply  that  the  mutation  effects 
will  be  beneficial  in  all  PTS  and  non-PTS  environments.  To  test  if  mutations  in  pykF  confer 
fitness  advantages  by  affecting  the  PTS  pathway,  we  measured  the  effect  of  one  of  the 
mutations  (pykF::\S150)  in  five  PTS  and  six  non-PTS  resources  in  the  genetic  background 
where  it  occurred  and  was  fixed  (Figure  14).  We  found  the  pykF::\S1 50  mutation  has  a 
pleiotropic  effect  across  environments:  it  is  neutral  in  most  of  the  PTS  environments  (mean  ± 
SEM:  0.981  ±  0.015),  but  has  different  deleterious  effect  in  most  of  the  non-PTS 
environments  (mean  ±  SEM:  0.877  ±  0.022).  A  two-way  ANOVA  with  environment  nested 
within  resource  type  (PTS  or  non-PTS)  found  a  significant  difference  between  PTS  and  non- 
PTS  resources  in  affecting  fitness  of  the  mutation  ( Fii32  =  302.9,  P  <  0.001 ).  This  result  is 
consistent  with  the  hypothesis  that  pykF  mutations  tended  to  confer  greater  fitness  in 
environments  requiring  PEP  for  uptake  than  in  alternative  substrates. 

Additional  work: 
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The  effect  of  a  mutation  may  depend  on  its  genetic  background  (epistasis),  the 
environment  (pleiotropy),  or  both.  However,  few  empirical  studies  have  investigated  the 
joint  influence  of  both  factors  in  the  same  experiment.  Those  that  have,  usually  focused  on 
randomly  generated  mutations,  which  may  behave  differently  than  beneficial  mutations 
selected  because  they  improve  organism  adaptation  in  a  specific  genetic  and  environmental 
context.  We  measured  the  fitness  effect  conferred  by  the  beneficial  pykF  deletion  mutation 
(A pykF)  in  23  natural  isolates  of  E.  coli  in  five  different  resource  environments  (Figure  17). 

We  found  that  the  genetic  background,  the  environment,  and  the  interaction  between  the 
genetic  background  and  the  environment,  all  have  significant  effects  on  the  fitness  effect  of 
the  mutation,  with  genetic  background  being  the  most  important  in  determining  fitness.  In 
agreement  with  the  diminishing  returns  epistasis  found  in  previous  experiments,  a  negative 
relationship  between  the  maximum  growth  rate  of  the  recipient  strains  and  the  fitness  effect  of 
Apy/cFwere  found  in  all  five  environments  (Figure  18),  although  only  in  one  environment  was 
this  relationship  significant.  Our  results  suggest  that  the  effect  of  beneficial  mutations  will  be 
difficult  to  predict,  depending  on  specific  interactions  with  genotype  and  environment. 
Nevertheless,  we  find  support  for  a  growing  body  of  research  predicting  a  qualitative 
relationship  between  mutation  effects  and  initial  fitness  of  a  progenitor  strain,  predicting  that 
the  contribution  of  a  mutation  will  be  smaller  when  added  to  fitter  progenitors.  (See  Fen 
Peng’s  PhD  Thesis  Chapter  4  for  details.) 
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Figure  14.  Fitness  effect  of 
pykF: :\S150  mutation  in  11 
resources.  Five  are  PTS  resources 
(red):  glucose,  fructose,  mannitol, 
mannose,  and  NAG;  six  are  non-PTS 
resources  (blue):  galactose,  glycerol, 
Resource  Type  lactose,  maltose,  melibiose,  and 
rhamnose.  The  neutral  effect  of 
pykF::\S150  in  glucose  seems 
contradictory  with  the  previous  6% 
benefit  shown  in  Figure  1,  because 
the  fitness  assays  here  were 
performed  in  96-well  plates  (see 
Materials  and  methods)  rather  than 
flasks. 
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Figure  15.  Box  and  whisker  plot 
indicating  PEP  concentration  of  otherwise 
isogenic  strains  expressing  the  indicated 
PykF  enzyme  ('Del1  is  the  deletion  allele). 

Strains  were  sampled  following  24  hours  of 
growth  in  the  evolution  environment.  Boxes 
indicate  first  and  third  quartiles,  the  central 
line  indicates  the  median,  and  whiskers 
indicate  95%  confidence  intervals  of  the 
median.  Symbols  at  top  of  panel  reflect 
results  of  Dunnett's  test  comparing  strains 
with  evolved  enzymes  against  the  WT 
enzyme  strain:  "***"  <0.001 ,  M*M  0.05, M  n  0.1 . 
Results  are  from  six  independent  replicate 
measurements  of  each  strain. 
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Figure  16.  Growth  and  fitness 
measurements  of  the  ancestral  strain 
and  derivatives  encoding  one  of  four 
different  pykF  alleles  over  24  hours. 

slipper  panel  shows  the  growth 
measurement  for  ancestral  strains  with  WT, 
Del,  and  A301S  allele,  and  error  bars 
indicate  95%  confidence  intervals  of 
independent  CFU  estimates  (n=6).  Lower 
panel  shows  the  fitness  effects  of  pykF 
alleles  at  different  time  points.  Symbols  and 
error  bars  represent  the  mean  and  95% 
confidence  intervals  of  four  independent 
measurements. 


Figure  17.  The  Apy/cF  mutation  was  added  to  the 
23  genetic  backgrounds  and  their  fitness 
measured  relative  to  the  original  progenitor  strain 
in  each  of  five  PTS  resource  environments. 

Hollow  symbols  represent  the  mean  fitness  of  the 
ApykF  mutation  in  each  strain  based  on  three  or  four 
independent  estimates.  Red  hollows  indicate  the 
fitness  of  ApykF  in  the  genetic  background  on  which 
it  arose  and  fixed.  Blue  squares  represent  the  grand 
mean  fitness  in  each  environment. 
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Figure  18.  Relationships  between 
maximum  growth  rate  and  relative 
fitness  of  adding  the  ApykF 
mutation  to  the  23  genotypes  in  five 
PTS  resource  environments. 

Hollows  represent  the  mean  of 
maximum  growth  rate  of  the  23 
genotypes  (3  <  n  <  32),  and  the  mean 
of  fitness  of  the  ApykF  mutation  in 
each  genotype  (n  =  3  or  4).  Lines 
show  linear  regression  with  P  value 
and  R2  calculated  in  a  Pearson 
correlation  shown  for  each  graph. 
Spearman  non-parametric  correlations 
were  also  performed,  and  they  did  not 
affect  conclusions  of  significance. 
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Impact 

Our  program  provides  a  uniquely  detailed  functional  understanding  of  how  evolution  by 
natural  selection  occurs  at  the  molecular  level.  Many  studies  have  described  the  genetic  and 
physiological  changes  that  occur  during  evolution.  Our  work  is  the  first  to  link  a  quantified 
fitness  advantage  to  specific  adaptive  mutations  via  changes  in  enzyme  structure  (dynamics), 
function  and  metabolic  consequences  in  the  context  of  the  E.  coli  long-term  evolution 
experiment. 

The  data  sets  generated  by  our  work  can  be  analyzed  in  the  context  of  a  wide  range  of  fields 
(e.g.  physiology,  biochemistry  and  bioinformatics).  Such  multidisciplinary  analyses  generate 
a  level  of  empirical  detail  that  will  be  invaluable  to  the  goal  of  developing  a  mechanistic 
understanding  and  prediction  of  the  process  of  adaptation. 

Independently,  this  work  also  provides  unique  insights  into  the  mechanism  of  pyruvate  kinase 
function.  In  particular,  we  have  been  able  to  delineate  how  FBP  binding  alters  the  structure 
and  dynamics  of  pyruvate  kinase  to  activate  the  enzyme.  These  data  has  been  lacking  for 
this  enzyme  until  now. 


Future  research 

We  envisage  two  future  programs  of  work: 

1)  We  have  demonstrated  that  the  adaptive  changes  in  pyruvate  kinase  are  mediated  through 
changes  in  protein  dynamics.  Is  this  a  general  phenomena?  The  gene  spoT  is  also  found  to 
contain  adaptive  mutations  in  most  of  the  12  populations.  Upon  starvation,  bacteria  slow 
down  the  physiological  processes  for  cell  growth  and  speed  up  the  processes  to  overcome 
nutrient  deficiency.  This  stringent  response  is  regulated  by  the  RelA  and  SpoT  enzymes, 
which  catalyze  the  synthesis  and  degradation  of  a  signaling  alarmone,  ppGpp.  Interestingly, 
the  SpoT  enzyme  of  E.  coli  comprises  both  a  ppGppase  and  GDP  synthetase,  thus  regulating 
the  ppGpp  levels  in  the  cell.  To  test  whether  the  adaptive  mutations  in  SpoT  also  alter  protein 
dynamics,  we  first  require  a  thorough  structural  and  biochemical  understanding  of  the  SpoT 
enzyme,  which,  unlike  pyruvate  kinase,  is  currently  lacking.  We  have  initiated  such  studies  of 
the  SpoT  enzyme  in  an  effort  to  understand  how  this  enzyme  functions  and  why  it  is  also 
target  for  adaptive  evolution. 

2)  During  the  course  of  our  program,  we  became  curious  about  the  nature  of  the  two  pyruvate 
kinase  deletion  mutations  that  are  found  in  the  long-term  E.  coli  evolution  experiment. 
Interestingly,  the  gene  deletions  result  in  a  truncated  (a/p)e-barrel  protein  -  the  fold  of  the 
catalytic  domain  -  such  that  the  resulting  protein  constitutes  a  half  and  quarter  barrel.  To 
date,  half  (a/(3)8-barrel  folds  have  only  been  engineered  and  have  never  been  observed  in 
nature  as  an  adaptive  change.  Moreover,  the  quarter  barrel  has  been  hypothesized,  but 
never  successfully  observed  (engineered  or  otherwise).  We  have  demonstrated  that  both 
truncations  can  be  overexpressed  and  are  soluble!  We  are  currently  working  to  purify  these 
proteins.  This  work  has  the  potential  to  significantly  add  to  the  story  of  how  (a/|3)8-barrels 
evolved. 
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